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Abstract The study concerns the question of how the shape sta¬ 
bility features of laminated columns of Norway spruce can be 
improved in terms of twist through optimal orientation of the in¬ 
dividual laminates. Both experimental testing and numerical si¬ 
mulations were employed for evaluating twist stability. In ah the 
columns studied, deformations were measured experimentally at 
different moisture content levels. A number of columns were 
also selected for numerical analysis in order to obtain a more 
thorough understanding of the twist behavior involved, their geo¬ 
metries and material properties of interest being determined. The 
experimental results showed the twist stability of the columns to 
be highly dependant upon the internal orientation of the indivi¬ 
dual laminates. It was also found that high quality columns in 
terms of shape stability could be manufactured, even when the 
center-core material has a strong twist tendency. The numerical 
simulations performed were in close agreement with the experi¬ 
mental results. 

Experimentelle und rechnerische Untersuchung 
der Formstabilitat von Kreuzholzbalken 

Zusammenfassung Diese Studie befasst sich mit der Frage, wie 
die Verdrehung von Kreuzholzbalken aus Fichte durch optimale 
Anordnung der zu verklebenden Teile reduziert und dadurch die 
Formstabilitat der Balken verbessert werden kann. Die Form- 
stabilitat der Balken wurde sowohl experimentell als auch mit 
Simulationsrechnungen untersucht. Bei alien untersuchten Bal¬ 
ken wurden Verformungen bei unterschiedlichen Holzfeuchten 
experimentell gemessen. Um das Verdrehungsverhalten besser 
zu verstehen, wurde dieses bei einigen Balken rechnerisch unter¬ 
sucht. Dazu wurde deren Aufbau und die relevanten Materialei- 
genschaften bestimmt. 
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Die Versuchsergebnisse belegten die groBe Abhangigkeit der 
Formstabilitat von der Anordnung der Einzelteile. AuBerdem 
zeigte sich, dass hochwertige formstabile Balken auch aus mark- 
nahem Holz mit starkem Drehwuchs hergestellt werden konnen. 
Die rechnerische Simulation stimmte mit den Versuchsergeb- 
nissen sehr gut iiberein. 


1 Introduction 

Distortions in sawn timber often come about through climatic va¬ 
riations or through growth stresses present in the log when it is 
sawn into planks or boards. These distortions are a serious pro¬ 
blem, both for the sawmill industry and for the building industry, 
since users often experience serious difficulties in handling dis¬ 
torted boards. It is not satisfactory for the user to get straight 
boards from the sawmill if these later become markedly distor¬ 
ted. Customers thus need timber and timber products that show 
a high degree of shape stability when exposed to climatic chan¬ 
ges, particularly since timber structures are normally subjected to 
considerable climatic variations frequently. The moisture-related 
deformations that occur are often divided into four types: twist, 
cup, spring and bow. Twist is the type of distortion considered 
to be the most severe. It represents the main cause of the re¬ 
jection of studs at building sites in Sweden (Johansson et al. 
1994, Perstorper 1994, Forsberg 1999, Johansson 2002). Du¬ 
ring the last decade, the use of engineered wood products in 
the building industry has increased. Most studies of shape stabi¬ 
lity have been concerned, however, with individual boards (e.g. 
Forsberg 1999, Woxblom 1999, Ormarsson et al. 1999, Serrano 
and Cassens 2000, Johansson 2002, Cooper and Maun 2003, 
Eriksson et al. 2004). The features of many composite wood pro¬ 
ducts are clearly in need of investigation. 

For laminated boards, it has been shown experimentally that 
gluing them in an appropriate way can substantially decrease the 
distortions that occur, particularly twist (see Serrano and Cassens 
2000, Cooper and Maun 2003, Eriksson et al. 2004). Studies that 
have been carried out have involved varying the orientation of the 
component parts, their dimensions and the tree species utilized, 
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such as Sitka spruce, yellow poplar and Norway spruce. Numeri¬ 
cal simulations have been performed earlier for investigating the 
influence of the material properties, cross-sectional size and an¬ 
nual ring orientation on the shape stability of individual boards; 
(see for example Ormarsson 1999, Ormarsson et al. 1999, Or- 
marsson et al. 2000). The model used in these simulations has 
been validated by studies of separate boards, (see Dahlblom et al. 
2000, Dahlblom et al. 2001), close agreement with experimen¬ 
tal observations being obtained. The moisture-induced distortion 
studies of laminated boards presented in Ormarsson (1999) and 
Ormarsson et al. (2001) concern numerical simulations, experi¬ 
mental support for the findings obtained numerically being pre¬ 
sented recently in Eriksson (2004) and in Eriksson et al. (2004). 
It was concluded that if the laminates are oriented in an optimal 
way, boards of excellent shape stability could be produced. Other 
timber products, such as columns and doors, have likewise been 
studied numerically in (Ormarsson 1999, Ormarsson et al. 2001, 
Ormarsson et al. 2002). The results obtained are of particular in¬ 
terest in the present investigation and will be taken as a starting 
point in examining the twist stability of columns made of four 
laminates. 

Some of the findings reported in Ormarsson (1999) are sum¬ 
marized in Table 1, which shows the degree and direction of twist 
that developed in columns of four types with the dimensions of 
3 x 0.1 x 0.1 m 3 , when dried from about 27% moisture content 
(MC) to 12% MC. The results are given in degrees of twist per 
running meter of the column (°/m). A square 35 x 35 mm 2 hole 


Table 1 Computed twist of various columns dealt with in Ormarsson (1999) 
Tabelle I Errechnete Verdrehung verschiedener Balken. Aus Ormarsson 
(1999) 

Name Cross- Twist A column exhibiting positive twist 

section (°/m) 



was made at the center of the cross-section in each column. The 
material data used throughout the simulations was that one typi¬ 
cal for Norway spruce under normal growth conditions. Case C4 
shown at the bottom has the largest twist of all the columns pre¬ 
sented in Ormarsson (1999), but is not part of this investigation, 
being included for comparison purposes only. It shows almost 
twice as much twist as Cl. Note that Cl and C4 are of opposite 
twist direction. An example of positive twist is shown in the dia¬ 
gram at the right in Table 1. Moving the pith in theoretical terms 
from its being in the corners in C1, to its being at the center of the 
product in C4 results in a —7.3 °/m change in twist. This large 
difference implies that products of this type can be very sensitive 
to changes in the pith position. Note too, that the cases C2 and 
C3, although consisting of the same material, show very good 
twist stability as compared with that of the Cl and C4-types. 


2 Material and method 

On the basis of previous findings, three different types of co¬ 
lumns were manufactured, the orientation of the laminates used 
in columns C1-C3 being employed, see Table 1, although the 
columns studied in Ormarsson (1999) had larger holes and so¬ 
mewhat larger cross-sections than those studied here. Each of the 
columns consisted of spruce timber (Norway spruce) stemming 
from the central part of Sweden. To investigate the twist stabi¬ 
lity of the manufactured products, both experimental testing and 
numerical simulations were employed. 

2.1 Geometry and sawing pattern 

Each of the products manufactured was cut to 2.4 m in length 
and was planed to 93 x 93 mm 2 in cross-section, an 18 x 18 mm 2 
square hole being made through the center of it. Each product 
stemmed from two center boards with split pith, as shown in 
Fig. 1. The sawn boards were cut through the pith once more to 
produce the four laminates that made up each column. On the 
basis of observations reported in Forsberg (1999) and Johansson 
(2002), it was concluded that boards sawn close to the pith deve¬ 
lop a roughly equal degree of twist when taken from the top end 
of a log as when taken from the butt end. Accordingly, no consi¬ 
deration was taken of a board’s longitudinal position in the tree. 
Note that C3 can be obtained from C2 by rotating and shifting 
the two lamellas on the right side. Moreover, columns of the Cl 
type, which are quite common on the market, can be manufac¬ 
tured from round timber, taking advantage of its wane when the 
square hole is made, thus saving material. 


Fig. 1 Sawing patterns and glued surfaces of the columns 
studied 

Abb. 1 Einschnitt und Anordnung der zu verklebenden 
Teile der untersuchten Balken 
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2.2 Industrial manufacturing of test specimens 

The experimental work was carried out in cooperation with 
a Swedish sawmill. The process of manufacturing the columns 
involved sawing and planing the timber to form the components 
and then gluing them together by use of a conventional gluing 
procedure. All the boards used in the study were piled in timber 
packages and were dried without top load in a low-temperature 
kiln until the MC had reached a level of about 15 to 18%. When 
the drying had been completed, the boards that had become most 
twisted were selected for manufacturing of the 20 Cl-type, 12 
C2-type and 8 C3-type products, respectively, that were to be 
tested. Finally, the products were wrapped in plastic to prevent 
further drying during their transport to the laboratory. 


2.3 Testing procedure 

After arriving at the laboratory, the columns were hung in a ver¬ 
tical position in a room of constant humidity and temperature, at 
85% RH and 20 °C, respectively. When the moisture content and 
the distortion had stabilized, the deformation of each column was 
measured. Displacements were recorded at the bottom and at one 
side, both in the middle and at the one end of the columns, the 
other end being fixed by a clamped support. A more thorough 
account of the frame setup and the equipment used to measure 
twist deformations can be found in Eriksson et al. (2004). The 
columns were stored then in another climate room with an RH of 
about 35% and a temperature of 20 °C for studying their shape 
stability. The experimental procedure just described was repea- 


Fig. 2 Changes in twist of the wooden co¬ 
lumns caused by a change in MC from 16.5% 
to 9.5%. Note that each pair of columns con¬ 
sists of parts from the same log. a Columns 
of the C2 type compared with those of the 
Cl type, b Columns of the C3 type compared 
with those of the Cl type 
Abb. 2 Veranderung der Verdrehung der 
Balken bei einer Holzfeuchteanderung von 
16,5% auf 9,5%. Hinweis: Jeweils beide 
Balken wurden aus dem gleichen Stammab- 
schnitt hergestellt. a Balken vom Typ C2 im 
Vergleich mit Balken vom Typ Cl. b Balken 
vom Typ C3 im Vergleich mit Balken vom 
Typ Cl 



1 2 3 4 5 6 7 8 9 10 11 12 


a Column number 



1 2 3 4 5 6 7 8 


b Column number 


Fig. 3 Twist deformations of the columns as 
measured at four different moisture levels (the 
bold marked columns were selected for nu¬ 
merical simulation studies, a Columns of the 
C2-type. b Columns of the Cl-type that stem 
from same log as the C2-type. c Columns of 
the C3-type. d Columns of the Cl-type that 
stem from same log as the C3-type 
Abb. 3 Verdrehung der Balken bei vier unter- 
schiedlichen Holzfeuchten (die fettgedruckten 
Nummern wurden fur Simulationsrechnungen 
ausgewahlt. a Balken vom Typ C2. b Balken 
vom Typ Cl und Typ C2 stammen aus dem 
gleichen Stammabschnitt. c Balken vom Typ 
C3. d Balken vom Typ Cl und Typ C3 stam¬ 
men aus dem gleichen Stammabschnitt 
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first moistening (16.5%) 
first drying (9.5%) 
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second drying (10.5%) 
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ted once. Results for the C3 products were compared both with 
those for the C1 and for the C2 products stemming from the same 
tree, see Fig. 1. 


3 Experimental results 

The major results of the experimental study are presented in the 
following. From the standpoint of shape stability, it is the change 
in twist between two different moisture states that is of primary 
interest. The absolute changes in twist when the columns were 
dried from 16.5% MC to 9.5% MC are shown in Fig. 2. Clearly, 
the C2 and C3 columns are much more twist-stable than the Cl 
columns. The mean change in twist was found to be 1.33 °/m for 
the Cl columns, and 0.44 °/m and 0.43 °/m, respectively, for the 
C2 and C3 columns. Thus, for the cases C2 and C3, the average 
change in twist was only about a third as great. Studying the ef¬ 
fect of repeated moistening and drying of the columns was also of 
interest. Looking more in detail on the results, the twist deformati¬ 
ons at equilibrium for four different MC levels are shown in Fig. 3, 
where positive twist is defined in the manner illustrated by the 
diagram in Table 1. As can be clearly seen, twist stability was ba¬ 
sically maintained during repeated moistening and drying. It can 
also be noted that for most of the C2 and C3 columns, drying led 
to a negative change in twist whereas for most of the Cl columns, 
the change in twist upon drying was positive. A few of the columns 
failed to follow this pattern. More specifically, the columns deno¬ 
ted as C2.1, C2.4, C2.8, C3.7 and C3.8 showed positive changes 
in twist upon drying whereas the corresponding Cl-type columns 
partly followed and partly failed to follow the general pattern de¬ 
scribed above. On the basis of the experimental results presented, 
two pairs of columns that deviated from the general twist pattern 
and two pairs of columns that followed it were selected for fur¬ 
ther study with the aid of numerical simulations. These pairs are 
marked with bold numbers in the diagrams. In addition, the de¬ 
formations measured were also used to obtain rough estimates of 
the bow at the middle of the columns. The average change in bow 
that moisture change led to was about 0.5 mm for the C2 and C3 
products and about 1 mm for the Cl products. 


Account was taken of the moisture and temperature dependency 
of the material parameters as well as the radial inhomogenity of 
the material properties and the spiral grain angle. Further theo¬ 
retical considerations of the model employed have been taken 
up (Ormarsson 1999, Ormarsson et al. 1999, Ormarsson et al. 
2000). All the columns were assumed to be stress-free at an in¬ 
itial MC of 14%. The products were moistened to a 16.5% MC 
and dried then to a 9.5% MC. Each of the columns was free to de¬ 
form during the variation in moisture, fixed boundary conditions 
only being introduced to prevent rigid body motions. 

4.1 Material data 

The most important material parameters for predicting twist sta¬ 
bility under varying climatic conditions are the cross-grain shrin¬ 
kage and the spiral grain angle. In addition, the pith position and 
the column geometry needed to be estimated. The test specimens 
were manufactured from the eight columns selected earlier. The 
origin of the specimens within each of the columns is illustra¬ 
ted in Fig. 4. The cross-sectional geometry and pith positions 
were assessed using a graded check-patterned slide, specimens 
for the shrinkage test being produced from clear-wood samp¬ 
les. To obtain estimates of variation in the spiral grain angle 
and in shrinkage, four and three samples, respectively, were ma¬ 
nufactured for determining these parameters. In cases in which 
the pith direction deviated from the longitudinal direction of the 
columns, this was compensated for during sample preparation. 
The values taken for the diffusion coefficients used in the si¬ 
mulations were those typical at 20 °C for Norway spruce, as 
given in the literature. The values chosen were 4 x 10“ 10 m 2 /s 
in the cross-grain directions and 13 x 10 10 m 2 /s in the longi¬ 
tudinal direction. Additional material data such as stiffness and 
mechano-sorption in Norway spruce were taken from Ormarsson 
(1999) as typical values. The coordinates used for defining the 
geometry of the columns were those representing the average of 
the coordinates as measured at the root end and at the top end of 
the various columns. Since each column consisted of four lamel- 
las, the pith position of each lamella had to be specified. The pith 
was assumed to vary in a linear manner between the coordinates 
at the two ends. 


4 Numerical simulations 

In order to obtain a more thorough understanding of the distor¬ 
tion behavior of the products that were studied, finite element 
simulations of the eight columns that had been selected were per¬ 
formed. The distortion model employed consisted of two parts, 
one for analyzing transient moisture flow and the other for ana¬ 
lyzing moisture-related deformations and stresses. In all the si¬ 
mulations, use was made of 30 x 16 x 16 elements. The sources 
of input to the moisture model were the initial moisture state, 
the corresponding equilibrium moisture content of the surroun¬ 
ding air, and the diffusion properties and orientation of the wood 
fibers. The distortions that were simulated were based on the 
moisture content histories obtained. The constitutive model con¬ 
sisted of elastic, moisture-induced and mechano-sorptive strains. 


Shrinkage properties. The samples were weighed continuously 
during moistening and drying. For each of the samples, moisture 
equilibrium was assumed to been reached when the change in 
weight was less than 0.001 g per 24 h. Shrinkage measurements 



Shrinkage 



Fig. 4 Origins within the columns of the material specimens of small size 
Abb. 4 Entnahmestellen der kleinen Proben aus dem Balken 
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were performed in a climate chamber under two different moi¬ 
sture conditions, 45% RH and 80% RH, corresponding to about 
10% MC and 14.5% MC, respectively, the dry density of the 
MC samples being determined later by use of the oven-dry me¬ 
thod (drying the sample for 24 h at 103 °C). The samples were 
200 x 10 x 10 mm 3 in size. A displacement gage was employed 
for measurements in the radial and tangential directions when 
the moisture in the samples had reached a state of equilibrium. 
The results are shown in Fig. 5, the averages obtained corre¬ 
sponding to typical values for Norway spruce as given in the 
literature, although the maximum and the minimum shrinkages 
that occur differ by a factor of about three in the radial direction. 
A high radial value appears to be accompanied by a high degree 
of tangential shrinkage and a low radial value by a low tangential 
shrinkage coefficient. On the average, the tangential shrinkage 
was about twice the radial. In Fig. 5, shrinkage parameters that 
pertain to the same column are denoted by a common symbol. 
The average of the three measurements obtained in each case 
is assumed to be representative for the column in question, the 
influence on the shrinkage parameters of the distance of the sam¬ 
ple from the pith being neglected. The means of the cross-grain 
shrinkage parameters as measured are listed in Table 2. 

Spiral grain. The spiral grain angle was determined using 
a slope-of-grain indicator that included a scriber and a protractor 
setup, the angle of the scriber groove being measured. The ave¬ 
rage of two spiral grain measures was obtained for each of the 
four laminates. To estimate the radial variation in the spiral grain 
angle, a few annual rings were peeled off after each reading. The 
measurement results are shown in Fig. 6. The average values ob¬ 
tained are similar to those given in literature. The coefficient of 
determination was found to be 0.12 for the linear variation of 
the spiral grain angle as seen in relation to the distance from the 
pith. This is comparable to what was found in a more extensive 
spiral grain study concerning Norway spruce, (see Perstorper et 
al. 2001). Although it may appear that the spiral grain angle is 
poorly correlated to the distance from pith when all the samples 
are examined at the same time, when each column is studied se¬ 


parately, the assumption of a linear relation appears justified, see 
Fig. 7. Note that the spiral grain angles shown in Fig. 7 vary rat¬ 
her much if one bears in mind the fact that all the material stems 
from the same tree. On the basis of the experimental observati¬ 
ons, the function for the spiral grain angle 6 (°) was taken to be 
linearly dependent upon the distance from pith r (mm), such that 

0 = A + Bxr (1) 

The coefficients A and B were determined by a least square fit of 
the measured data to Eq. 1, see Table 2. 

4.2 Simulation results 

The numerically simulated twist results together with the mea¬ 
sured ones are shown in Fig. 8. In most cases, the numerical 
predictions agree well with the experimental findings. The be¬ 
havior of the twist stable column C2.6 and the corresponding 
twist-prone one column Cl.6 is closely predicted by the nume¬ 
rical results. Simulation showed that the twist stability of C3.7 
and C1.7 was approximately equal and the experimentally mea¬ 
sured twists were in reasonable accordance. In the simulations, 
C3.3 showed only a small change in twist upon drying in con¬ 
trast to Cl.3. The experimental results support the behavior of 
both types, the deviation from the predicted twist value being 

Table 2 Coefficients for the spiral grain function and the cross-grain shrin¬ 
kage parameters as measured for the eight columns selected for numerical 
simulations 

Tabelle 2 Koeffizienten der Gleichung fur den Faserwinkel sowie gemessene 
radiale und tangentiale SchwindmaBe der acht fur die Simulationsrechnung 
ausgewahlten Balken 



C2.1 

Cl.l 

C2.6 

Cl.6 

C3.3 

Cl.3 

C3.7 

Cl.7 

A 

3.79 

4.02 

5.43 

4.22 

7.30 

6.21 

1.38 

2.28 

B 

-0.131 

-0.056 

-0.043 

0.006 

-0.086 

-0.067 

-0.034 

-0.058 

a r 

0.106 

0.199 

0.164 

0.230 

0.176 

0.171 

0.186 

0.214 

a, 

0.290 

0.308 

0.295 

0.330 

0.384 

0.314 

0.339 

0.365 



0.40 

0.36 

0.32 

0.28 

0.24 

0.08 0,12 0,16 0,20 0.24 0.28 





Distance from pith (mm) 


Fig. 5 The radial versus tangential shrinkage parameters. Each symbol cor¬ 
responds to one particular column and the solid lines denote the overall 
average 

Abb. 5 Zusammenhang zwischen radialen und tangentialen SchwindmaBen. 
Jedes Symbol entspricht einem bestimmten Balken. Die durchgezogenen 
Linien geben die Gesamtmittelwerte an 


Fig. 6 Spiral grain angle versus distance from the pith for eight tested co¬ 
lumns. The solid lines denote the average linear relation and the cross marks 
the measured value of the spiral grain angle 

Abb. 6 Zusammenhang zwischen Faserwinkel und Abstand vom Mark von 
acht untersuchten Balken. Die durchgezogenen Linien geben die mittlere 
Abhangigkeit an. Die Kreuze den einzelnen Messwert 
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Fig. 7 Spiral grain angle versus distance from the 
pith for column C2.1 and the corresponding Cl- 
type column 

Abb. 7 Zusammenhang zwischen Faserwinkel 
und Abstand vom Mark des C2.1 Balkens und 
des entsprechenden Balkens vom Typ Cl 




Distance from pith (mm) 


greater in the latter case. The column C2.1, which showed an 
unexpectedly large twist during the experiments, was found to 
behave in basically the same way in the simulations. The beha¬ 
vior of the corresponding Cl.l column, however, was not pre¬ 
dicted in a satisfactory way, showing numerically a change in 
twist of opposite sign to that obtained in the experiments. The 
spiral grain function was found to be quite unlike that for C2.1 
and for the corresponding C1 column, although the lamellas stem 
from same tree, see Fig. 7. It could have been expected that the 
spiral grain function for the various products would have been si¬ 
milar, yet spiral grain measurements are generally difficult and 
may be sensitive to local anomalies in fiber directions, such as 
through closeness to a knot. The large differences observed and 
the close agreement with the results of the experiments motiva¬ 
tes use of the function obtained for C2.1 in computation of the 
Cl product. The twist became -0.28 °/m instead of 0.86 °/m, the 
former being in closer agreement with experimental findings, see 
the striped bar in Fig. 8. In a manner similar to what was found 
for C2.1, Cl.3 showed a somewhat larger twist in the simulations 
through use of the spiral grain function obtained for C3.3. 


5 Conclusions and discussion 

On the basis of results presented it can be concluded that for a gi¬ 
ven change in moisture content the C2 and C3 columns were 
more twist-stable than the corresponding Cl columns. The co¬ 
lumns were shown to differ in the direction of twist obtained 
upon drying despite their consisting of similar material. During 
repeated moistening and drying, all the columns maintained their 
shape stability in terms of twist, there being no evidence of pro¬ 
ducts tending to deviate from the twist pattern found generally. 
None of the products studied showed any substantial moisture- 
induced spring or bow. In most cases, results for simulated twist 
agreed well with the corresponding experimental results. The 
C2.1 and the corresponding Cl product deviated from the gene¬ 
ral twist pattern to some extent. Numerical simulation showed 
that this could be explained, however, by the pith being displa¬ 
ced from its normal position as well as by the large slope of spiral 
grain angle. Deviations in the twist pattern of the C3.7 and Cl.7 
columns were found to mainly be due to the spiral grain function 
involved since the pith was located close to its intended position. 
The same can be thought to apply to C3.8, the spiral grain func¬ 
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Fig. 8 Comparison of numerical and experimental results for twist for a se¬ 
lected set of columns following a change in MC from 16.5 to 9.5% 

Abb. 8 Vergleich der rechnerisch und experimentell bestimmten Ergebnisse 
der Verdrehung ausgewahlter Balken bezogen auf eine Holzfeuchteanderung 
von 16, 5% auf 9, 5% 


tion of which was found to deviate markedly from the average, 
although this product was not studied further. 

When the correct shrinkage parameters, the spiral grain an¬ 
gle and the pith position were taken into account, the simulation 
tool was found capable to predicting well the twist which ac¬ 
tually occurred in the individual columns and whether or not 
it would obey the general pattern of twist. It appears therefore, 
that products sensitive to deviations in pith location should be 
analyzed individually with the correct pith position being ta¬ 
ken into account. In such an analysis, one needs to clarify how 
the properties of the material, its internal structure and the en¬ 
vironmental conditions affect the moisture-related deformation 
process. It can be concluded that the technique developed here 
for orientating the various component laminates relative to each 
other can be used to reduce twist deformation markedly in other 
types of products as well since columns with good shape stability 
are obtained, even from center-core timber with a strong twist 
tendency. 
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